New scenario for high-T c cuprates : electronic topological transition as a motor for 

anomalies in the underdoped regime. 
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We have discovered a new nontrivial aspect of electronic topological transition (ETT) in a 2D free 
fermion system on a square lattice. The corresponding exotic quantum critical point, 5=S C , T=0, 
(n=l — S is an electron concentration) is at the origin of anomalous behaviour in the interacting 
system on one side of ETT, 8 < 8 C . The most important is an appearance of the line of characteristic 
temperatures, T*(S) oc S c — 5. Application of the theory to high-T c cuprates reveals a striking 
similarity to the observed experimentally behaviour in the underdoped regime (NMR and ARPES). 
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This is a particularly exciting time for high-T c . The 
experimental knowledge converges. Almost all experi- 
ments, NMR Kg], ARPES @, infrared conductivity |] 
etc. provide an evidence for the existence of a characteris- 
tic energy scale T* (6) in the underdoped regime (S is hole 
doping). Below and around the line T*(S) the "normal" 
state (i.e. above T c ) has properties fundamentally incom- 
patible with the present understanding of metal physics. 
The field has reached the point when a consistent theory 
is necessary to understand this exotic from theoretical 
point of view but quite well defined from experimental 
point metallic behaviour. The issue has a significance 
beyond the field of high-T c superconductivity - the funda- 
mental question arises : what kind of metallic behaviour 
is there, in addition to the well understood Fermi liquid? 

In the paper we propose our variant of the answer. We 
reexamine a free electron 2D system on a square lattice 
with hopping beyond nearest neighbors. We show that 
when varying the electron concentration defined as 1 — S, 
the system undergoes an electronic topological transition 
(ETT) H at a critical value 5=S C . The corresponding 
T=0 quantum critical point (QCP) combines several as- 
pects of criticality. The first standard one is related to 
singularities in thermodynamic properties, in density of 
states at w=0 (Van Hove singularity), to additional sin- 
gularity in the superconducting (SC) response function 
(RF) M. The second nontrivial aspect is: the same QCP 
is the end of the critical line T=0, 6 > 3 C each point 
6 of which is characterized by static Kohn singularity 
(KS) in polarizability of 2D free fermions. [What we 
mean as a static KS is a singularity at the wavevector 
connecting two points of Fermi surface (FS) with paral- 
lel tangents H], The two aspects of criticality are not 
related with each other. It is the latter aspect (never 
considered before) which as we will show is a motor for 
anomalous behaviour in the regime < 5 < S c of the 
system of noninteracting and interacting electrons (or of 
any fermion-like quasiparticles e.g. of those Q appearing 
in the t — t' — J model describing the strongly correlated 
CuOi plane responsible for main physics in the cuprates) . 
The found anomalies have a striking similarity to anoma- 
lies in the underdoped high-T c cuprates. The effect exists 
in all cases t' ^ or/and t" ^ ... except for special sets 



of the parameters corresponding to the perfect nesting in 
FS (including i' =£"=...— > 0) studied in many papers, see 
e.g. ||. For such sets the QCP loses the latter aspect of 
criticality and the anomalies disappear. 

A starting point is a 2D electron system on a square 
lattice with hopping beyond nearest neighbors 



£k = — 2i(cos k x + cos k y ) — At' cos k x cos k y 



(1) 



For any set of the parameters t, t', t" the dispersion 
law is characterized by two different saddle points (SP's) 
located at (± n, 0) and (0, ±7r) with the energy e s . When 
we vary the chemical potential fj, or the energy distance 
from the SP, Z = fi — e s , the topology of the FS changes 
when Z goes from Z > to Z < through the critical 
value Z — 0. In vicinities of SP's the dispersion law is : 



e(k) = e k — (i = -Z + ak^ — bkp, 



(2) 



where k is measured from (0, tt) (a = x, f3 = y) or from 
(ir, 0) (a = y, (3 — x). Explicit expressions for a and b 
depend on t, t',... We consider the general case : a ^ 0, 
b 7^ 0, a 7^ b. We choose a > b corresponding to t'/t < 0. 

The T = ETT has several characteristic aspects. The 
first (trivial) one is related to the local change of FS 
topology in the vicinity of SP. This leads to divergences 
in thermodynamic properties, in density of states at u> = 
etc. From this point of view the corresponding QCP is 
of a gaussian type with the dynamic exponent z = 2. 

The nontrivial aspect is related to mutual change in 
topology of FS in vicinities of two different SP's and re- 
veals itself when considering the electron polarizalibility 



n (ek) - n" (e q+k ) 



q+k — £k 



(3) 



We show that the latter has a square-root singularity at 
lu=0 and wavevector q=q m in a vicinity of Q=(7r, tt) for 
any Z on the semiaxis Z < 0: x°(q; 0) — x°(q m , 0) oc 
Vrlm q] for q < q m - ft is a static KS in the 2D elec- 
tron system. The locus of the wavevectors q m in the BZ 
is a closed curve around Q with |Q — q m |oc y/\ Z \. With 
decreasing | Z \ the closed curve shrinks and is reduced to 
the point q = Q at Z = where x O (q>0) diverges log- 
arithmically. The curve of the static KS's with q close 
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to Q does not reappear for Z > 0: X^IjQ) is peaked 
at q=Q in an intimate vicinity of ETT and it exhibits 
a wide plateau around q = Q for larger Z. To illustrate 
this we show in Fig.l the q dependence of x°(qi0) cal- 
culated based on (|J) and ([!]). [We use the model with 
only t'^= being a generic model for the family: a=/= 0, 
0, b] The discussed above curve is the closest to 
q = Q curve of singularities in Fig. la. In the plots one 
sees only a quarter of the picture around q = Q; to see 
the closed curve around (it, tt) one has to consider the 
extended BZ. [Few other curves of KS's seen in Fig.l are 
not sensitive to ETT, we discuss them elsewhere.] 




FIG. 1. q dependence of X°(ci>0) through the Brillouin 
zone for (a)Z <0, (b)Z >0, (c)Z=0. Q x =q x /ir, Q y =q v /^. 
The point q = Q corresponds to the left corner, (t 1 /t — —0.3) 

As a result, the point Z = 0, T = turns out to be 
the end point of the critical line Z < 0, T = 0. 

Paradoxically, the absence of the discussed curve of 
static KS's for Z > leads to anomalous behaviour 
of the system on this side of QCP. To see this let's 
calculate u> dependencies of i?ex°(q, u>), /mx°(q, oj) and 
C(ui)=Imx (q, o-i)/u> for the characteristic for this regime 
wavevector q=Q. Results are shown in Fig. 2a. One can 
see that all functions are singular at some energy uj c . 

Calculations with the hyperbolic spectrum (|2|) give the 
following expression: 

Re X °(Q,uj) = Rex°{Q,u c ) - $(u/u c ,b/a)/t with 



F(x,y) 




7i(2/)(l-x 2 ), 

i2(y)Vx- l, 



x-l< 
< x- 1< 1 



(4) 



(71,2(2/) *C 1). The new energy scale which appears 
and corresponds to the singularities in Fig. 2a is given by 



lu c = Z(l + b/a). 
The singularities at oj — uj c are dynamic 2D KS's. 
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FIG. 2. Rex°(Q,u) (solid), 7mx°(Q,w) (dot-dashed) and 
C(tu)—Imx (Q, lu)/uj (dashed) in the regime Z > (a) as a 
function of uo for T = 0, (b) as a function of T for ui — ► 
(7mx°(q, 0) = 0). Here t'/i = -0.3, Z/t = 0.21. 

The dynamic KS's at T — are transforming into 
static Kohn anomalies at finite temperatures which again 
scale with Z, see Fig. 2b. When comparing with Fig. 2a 
one can see that the behaviour is similar being smoothed 
by the effect of finite T (both behaviour are very dif- 
ferent from the standard one for noninteracting system). 
The important difference is that the characteristic tem- 
peratures of the Kohn anomalies for i?ex°(Q,0) and for 
lim^^o imx°(Q, uj)/lo being both scaled with Z 



T* Re = AZ, 



A < B 



are different on the contrary to the characteristic energy 
of the KS's at T=0 (that is an usual effect of finite T). 

Another remarkable signature of asymmetry in Z 
is the following. Taken for the characteristic for each 
regime wavevector, q=q m for Z< and q = Q for Z > 0, 
X°(q, 0) decreases rapidly with \Z\ for Z< (ordinary 
behaviour) while for Z> it remains practically con- 
stant (and quite high) for Z not too small. Moreover 
for finite T, x°(Q,0) has a maximum at Z=Z*(T)> 
(that is very unusual). As a result of the described T 
and Z dependencies of x°(Qj 0) in the regime Z > 0, the 
lines X°(Q) 0)=const have an unusual form in the T — Z 
plane: they develop rather around the "critical" lines 
T* Re (Z) and T* m (Z) than around the QCP, T=0, Z=0. 

On the contrary, a behaviour of SC RF (in both cases 
isotropic s-wave or d-wave symmetry) is symmetrical in 
Z being related to the first aspect of ETT. For the same 
reason the SC RF decreases quite rapidly with increasing 
a distance from QCP, i.e. with increasing T and \Z\. 

Above we considered a system of noninteracting elec- 
trons. In fact the same picture takes place for any system 
of fermion or fermion-like quasiparticles which dispersion 
law is determined by the topology of 2D square lattice 
and has a form ([!]) . In || where we discuss some problems 
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of strongly correlated systems we show that such quasi- 
particles (with spin and charge) do exist in the t — t' — J 
model describing the strongly correlated CuOi plane. On 
the other hand, the shape of FS observed by ARPES does 
imply the existence of nnn hopping t' ^ 0, so that the 
condition of the asymmetry a ^ b necessary for the ex- 
istence of the discussed ETT is fulfilled. Moreover this 
shape implies t'/t < 0, the case for which the critical 
doping S c is positive. Below we will pass from the energy 
distance from ETT, Z, to the doping distance, S c — 5, 
using large FS condition: 1 — S=2J2] SL nF (^k) @- 
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FIG. 3. Phase diagram with the lines of SDW and d-wave 
SC instabilities and the lines T^ e (S), Tj m (S) (t'/t = -0.3, 
t/J = 1.9). We consider only the metallic part of phase dia- 
gram (for discussion about a passage from AF localized-spin 
state at low doping to the metallic state with large FS for 
intermediate doping see [8]). 

Let's consider now the system in the presence of in- 
teraction. A quite trivial consequence of the ETT is a 
developing of density wave (DW) and SC instabilities 
around the point S=S C , T=0. [The effects are related 
to the logarithmic divergence of x°(QiO) an d lnZlnT 
divergence of the SC RF as T -> 0, Z -> 0.] Non- 
trivial consequences concerning the DW degrees of free- 
dom and related to the Kohn singularity aspect of ETT 
are: (i) strong asymmetry between regimes d < S c and 
S > 5 C and (ii) very long (in doping and temperature) 
memory about DW instability in the disordered state on 
one side of ETT, 5 < 5 C . To see this let's consider the 
electron-hole RF which in the simplest RPA approxima- 
tion is given by x(q,^) = X°(q,w)/(1 + V r q x°(q,w)). 
In the case of interaction Vq in a triplet (singlet) chan- 
nel the instability and normal state fluctuations are of 
SDW (CDW) type. We will consider the former inter- 
action: Vq—Jq—2J(cosq x + cosq y ) (J > 0) as strongly 
supported by neutron scattering and NMR experiments 
for the cuprates and on the other hand as an interaction 
between the discussed above quasiparticles in the t—t' — J 
model (|]. For such interaction both instabilities d-wave 
SC (see details in || ) and SDW take place around QCP. 
Due to the symmetry of SC RF in Z, T sc (5) is symmet- 
rical on two sides of S c with a maximum at <5 = 5 C , see 
Fig. 3. Therefore the regimes S < 5 C and S > S c can be 
considered as underdoped and overdoped, respectively. 
On the contrary, the line of SDW instability, Tsdw($)> 



given by x°(q, 0) = -1/ J q (q = Q for 6 < S c and q = q m 
for 6 > S c ) has an anomalous form in the regime 5 < 6 C : 
it develops rather around the lines T Re (S) and Tj m (S) 
than around QCP (see Fig. 3) reproducing the form of 
lines X°(Q>0) = const discussed above. 

When at certain doping, 8=8 sdw, the ordered SDW 
solution disappears, it is the disordered metallic state 
which keeps this type of behaviour: the regime T Re (8) < 
T < Tj m (5) (II) turns out to be a regime of a minimum 
disorder and the regime T < T Re (8) (I) a regime of a 
reentrant in temperature quantum SDW liquid. 
Indeed, two most important parameters characterizing 
SDW liquid: k 2 = 1 - |Jq|x°(Q,0) describing a "prox- 
imity" to the SDW instability and Tq = k 2 /C(0) de- 
scribing a relaxation energy behave in a reentrant way in 
increasing T: n 2 decreases (slightly) with T until T Re (8) 
and Tq decreases (strongly) until T* Re {8) < T£ < Tf m (S) 
as if the system would move towards an ordered phase. 
However, it does not reach it, the reentrancy stops and 
the system passes to the regime II of a minimum disor- 
der above which a standard disordered state behaviour 
is restored (regime III). On the other hand, the quantum 
SDW liquid state in the regime I is practically frozen in 
doping due to the very weak dependence of k 2 on dop- 
ing. As a result the disordered metal state in the regime 
5 < S c keeps a strong memory about the ordered SDW 
phase (and therefore develops strong critical SDW fluc- 
tuations) very far in doping and in temperature. On the 
contrary, in the regime S > S c the memory about SDW 
instability and the corresponding fluctuations disappear 
rapidly due to the sharp decrease of x° (q ; 0) with in- 
creasing 8 — S c and T. The same is valid in both regimes 
S > S c and S < 5 C for SC fluctuations due to the dis- 
cussed in the first part behaviour of SC RF as a function 
of T and \Z\. Therefore, although the SDW phase itself 
is energetically unfavorable with respect to the SC phase 
(except of the case of very high J/t), the metal state 
above T sc in the underdoped regime is a precursor of the 
SDW phase rather than of the SC phase. 
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FIG. 4. 1/TiT and 1/T 2G (a) calculated for S = 0.15 
(t'/t — —0.3, t/J — 1.9) (should be considered only above 
T sc ) and (b) taken from NMR for YBCOe.a [2]. 

The lines T Re (S) and Tj m (5) are basic lines for anoma- 
lies in the disordered metallic state. To demonstrate how 
the anomalies appear for different properties we consider 
some examples. In Fig. 4a we show calculated quasistatic 
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magnetic characteristics corresponding to the measured 
by NMR 1/TiT and 1/T 2G on copper as functions of T. 
A physical reason for a slight increase of 1 /T2G extending 
until w Tft e and a much stronger increase of \jT\T ex- 
tending until T ~ Tp is the discussed above reentrant be- 
haviour of k 2 and Tq with T. The theoretical behaviour 
is very close to that observed experimentally (Fig. 4b) and 
explains it in fact for the first time. 

In Fig. 5 we show an electron spectrum calculated for 
the ordered SDW phase (a) and for the disordered metal 
state (namely for the regime II) (b). For the ordered 
phase the spectrum is given by: £1,2 = (ca + e s)/2 ± 
y/(e A - es)/2) 2 + A 2 (e A (k) = e(k), e B (k) = e(k + Q)) 
with the gap A determined selfconsistently in the usual 
way. For the disordered state the " spectrum" is obtained 
from maxima of electron spectral functions strongly 
renormalizcd due to interaction with the described above 
SDW fluctuations. The characteristic form of the spec- 
trum in both cases is a result of a hybridization of two 
parts of the bare spectrum in the vicinity of two different 
SP's (0, 7r) and (tt, 0). The hybridization is static for the 
ordered SDW phase and is dynamic for the disordered 
state. [Details about the pseudogap opening in the dis- 
ordered state and its behaviour with T and S will be a 
subject of a separate paper.] The spectrum is in excel- 
lent agreement with ARPES data, see Fig. 5c (ARPES 
measures only the part corresponding to ui < 0). The ef- 
fect of splitting into two branches and of the pseudogap 
disappears quite rapidly in the regime 5 > 6 C due to the 
rapid weakening of SDW fluctuations. Due to the same 
reason it disappears roughly above Tf m (S). Both facts 
agree with experiments for the cuprates. 

SDW Disordered Experiment 
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FIG. 5. Electron spectrum e(k)/t along F — X symmetry 
lines, (a) in SDW phase (Z/t=0.d3 (5=0.25), T=0), (b) in 
the metallic state above T sc (Z/t=0.3 (5=0.1), T/t=0.15), 
(c) ARPES data [10] for underdoped Bi2Sr2CaCii20$+deita 
above T 3C . The dashed lines correspond to the bare spectrum, 
the thin line in (a) to the spectrum with the spectral weight 
less than 0.1. t' /t—— 0.3, t/J—1.8, wavevectors are taken in 
unit of 7r. 



as k 2 almost does not change with 5, the position of the 
peak does not as well. This agrees with INS data and 
explains (for the first time) the existence of the charac- 
teristic energy (~ 30 meV) above T sc for all 5, see e.g. the 
summarizing picture in Fig. 25 in JTT[ . As was emphasized 
before, strong SDW fluctuations disappear in the over- 
doped regime 5 > S c . In the underdoped regime they dis- 
appear (or strongly diminish) above Tf (5). Both facts 
are in a good agreement with INS. 

Summarizing, the simple picture arising from the ef- 
fect of ETT in 2D electron system on a square lattice 
gives a unified vision of normal state anomalies in the 
underdoped high-T c cuprates for both magnetic and elec- 
tronic properties. We succeed to explain the temper- 
ature anomalies in 1/T\T and I/T^g NMR characteris- 
tics, some crucial features of INS in the normal state, the 
disappearance of magnetic fluctuations in the overdoped 
regime, an opening of a pseudogap in electron spectrum, 
the shape of the latter in a vicinity of (0,7r), the disap- 
pearance of the pseudogap in the overdoped regime. All 
these are most nontrivial experimental results. Regard- 
ing that the theory does not use any external phenomcno- 
logical hypothesis and only two microscopical parameters 
t' /t and t/ J, the similarity between the theoretical results 
and experiments seems quite remarkable. We emphasize 
that the effect exists for any t'/t, t"/t,... except of two 
limit cases: (i) isotropic one a—b in eq.(2) (t'=t"=...=Q) 
and (ii) extreme anisotropic one a=0 or 6=0. Although 
ETT exists in both cases, the corresponding QCP's be- 
long to different classes of universality. For a=fe (nesting) 
the behaviour is symmetrical in Z, the anomalous regime 
discussed in the paper disappears. 
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Let's discuss a behaviour of 7mx(q,w), the character- 
istics measured by inelastic neutron scattering (INS) . As 
follows from the previous analysis, below Tj m it has a 
maximum at w=uo<xk 2 (being peaked at q=Q). So far 
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